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We analytically and numerically investigate the evolutions of the scalar perturbations through 
the cycles with nonsingular bounce. It is found that the amplitude of the curvature perturbation 
on large scale will be amplified cycle by cycle, and the isocurvature perturbations also obtain an 
amplification, but the rate of its amplification is slower than that of curvature perturbation, unless 
its coupling to the metric perturbation is not negligible. 



(N 

o 

(N 
C 



o 
cr 



O 



X 



I. INTRODUCTION 

In recent decades, inflationary cosmology has achieved 
a lot of successes, namely, it solved the homogeneity, 
isotropy, flatness problems and offered a noble mecha- 
nism of generating a nearly scale invariant spectrum of 
primordial fluctuations. However, as the inflationary cos- 
mology has singularity problem, it still deserves a further 
investigation on relevant issues. 

A cosmological cyclic scenario was proposed long time 
ago [l[ . In these years it has become popular 1 , which of- 
fered a new insight into the origin of observable universe. 
The inflationary model bases on the idea that space and 
time had a beginning when the universe had nearly in- 
finite temperature and density. In the cyclic model, the 
universe evolves cyclically in time, undergoing endless ex- 
pansions and contractions, cooling and heating, in which 
the density and temperature remain finite through these 
cycles. 

There have been many studies on cyclic models 

1,1, 1,1,0, 1,1, 0,[l3,[ll : [ll,p, in which it is 
generally assumed that the cyclic universe is homoge- 
neous cycle by cycle. However, when the perturbations 
are considered, the case mig ht be altered[Ii,[I3,[l3]. In 
general, the amplitude of curvature perturbation on large 
scale will increase in every contracting phase on large 
scale, however it will be nearly constant in the expanding 
phase. Thus the net result of one cycle is that the ampli- 
tude of curvature perturbation is amplified. This ampli- 
fication will possibly continue cycle by cycle. In certain 
sense, the amplification of curvature perturbation might 
eventually destroy the homogeneity of background. The 
effect of the increasing of other perturbations on cyclic 
universe has been also studied, e.g. anisotropic pressures 

m- 

Since the perturbation modes is continuously increas- 
ing through cycles, the reliability of some models of cyclic 
universe might be required to reevaluate. In the model 
building of a consistent cyclic universe, the perturbation 
should not only satisfy the requirement of the observable 
universe, but also is not expected to destroy the homo- 
geneity of background. Thus the study of the evolutions 
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FIG. 1: The upper panel is the evolution of scale factor in the 
model with increasing cycles. The lower panel is the evolu- 
tion of scale factor in the model with equal cycles. The time 
is scaled by k s which is an adjustable parameter to make the 
models be able to be consistent with the time scale of the 
observable universe. The insets are the evolution of u) corre- 
spond to every cyclic model around the bounce time. 



of cosmological perturbations in the model of cyclic uni- 
verse is significant. 

In this paper, we firstly recheck the evolution of cur- 
vature perturbation on large scale both analytically and 
numerically. It is proved that the curvature perturba- 
tion and the metric perturbation will be amplified cycle 
by cycle. Then we calculate the evolution of the pertur- 
bation of a light scalar field through cycles, which con- 
tributes the isocurvature perturbation. It is found that 
the isocurvature perturbations also obtain an amplifica- 
tion, but the rate of its amplification is slower than that 
of curvature perturbation, unless its coupling to the met- 
ric perturbation is not negligible. 

The outline of paper is as follows. In section II, we in- 
troduce the models of cyclic universe, which will be used 
in following sections, and the basic equations of perturba- 
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tions. In section III, we show the evolutions of the curva- 
ture perturbation on large scale through cycles, and also 
the metric perturbation, which will be used in section IV, 
in which we study the behavior of the isocurvature per- 
turbation on large scale without or with the effect of the 
metric perturbation. The final section is the discussion. 



II. THE SETUP AND THE BASIC EQUATION 
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FIG. 2: The evolution of background field (j>. The upper 
panel and the lower panel are these of the models in Figfj] 
corresponding with increasing cycles and equal cycles, respec- 
tively. 

In this section, we will introduce the models of cyclic 
universe, which will be used in following sections, and 
the basic equations of perturbations. We will regard the 
turnaround time as the beginning of a cycle, which is tL 
in the j th cycle. In each cycle the universe will sequen- 
tially experience the contraction, bounce, and expansion, 
and then arrive at the turnaround, which signals the end 
of a cycle. The beginning and the end of the contract- 
ing phase and the expanding phase in the j th cycle are 
jj fj fj 

l Cv u Ce> l Ei 



and vL , respectively, and the bounce is 



t B . Throughout this paper, "the model with equal cy- 
cles" means that the scale factor is periodically equal, 
while "the model with increasing cycles" means that the 
maximal value of the scale factor or the scale of cycle is 
cyclically increasing. 

The upper panel in FigQ]is the evolution of scale factor 
of the model with increasing cycles, e.g. see @ for details. 
The models with increasing cycles have been obtained in 
e -g-Hl)13|,[iI|> in which the entropy is increased, however, 
in our paper it results from a rapid expanding period after 
the bounce in each cycle, the corres ponding observable 
signals have been studied in e.g. [l9j.[20|.[2lj. The lower 
panel is that with equal cycles. We also give the evolution 



of background field 4> in Figj2] 

In Fig.l, the nonsingular bounce is implemented by 
applying a ghost field. In general, the introduction of 
ghost field is only the approximative simulation of a fun- 
damental theory below certain physical cutoff. In this 
sense, its appearance is only an artefact of this approx- 
imation. There have been lots of studies on how the 
perturbations pass through such a nonsingular bounce 
iH,m,[2i,[2l,[2i|- Recently, the nonsingular bounce 
has been obtained in higher derivative theories of gravity 
[13, or Galileon field [H|,|3(|, which areghost free, 



and others [3l|,[32|], and also earlier e.g.[33|.[34j|. 

We will investigate the evolutions of the scalar pertur- 
bations through cycles, and will not involve the details of 
model building. The models in FigQ] will serve the pur- 
pose of numerical simulations in the calculations of per- 
turbations. We shall work in longitudinal gauge in which 
if there is not the anisotropy only the metric perturbation 
<I> exist. The perturbation equations are 35], 36], [37] 
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where Vi denotes the derivative of the scalar field po- 
tential with respect to <f>. In principle, the modifica- 
tions implementing bounce are only reflected on the per- 
turbation modes at UV scale. While the perturbation 
modes we observed in today's cosmological observations 
are mainly in far infrared regime, thus this set of equa- 
tions are trustable. 

We, with this set of equations and the background evo- 
lutions in Fig.l, will analytically discuss and numerically 
simulate the evolutions of the curvature perturbation and 
the isocurvature perturbation on large scale in the follow- 
ing sections. 



III. THE EVOLUTION OF CURVATURE 
PERTURBATION 

We will recheck the evolution of the curvature per- 
turbation. The background field is <f>, and other fields 
contributing the isocurvature perturbation are neglected. 
Thus the curvature perturbation is ( = $ + ^S(/>. The 

equation of curvature perturbation in the momentum 
space is [38] , [11] , [Zo[ 



k 2 - — )u k = 0, 



(4) 
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where Uk — zQk and z — cup/H, the prime denotes the 
derivative for the conformal time rj. 

When the perturbations are deep inside the Hubble 
scale k 2 3> z"/z, the solution of perturbation is Uk ~ 
— h=e~ lkr * . While the perturbations are outside of Hubble 

scale k 2 <C z" /z, Eq.flU has general solution Uk, which 
gives e.g.jioj. 



Ck ~ C( 



C2 



dq 

72 



(5) 



where C^i and C^2 are constant only dependent of k. 
The perturbation spectrum will be determined by which 
one of both terms is dominated. 

In general, for the contraction phase 



a 



(6) 



where t < ts and n is constant. Thus during the contrac- 
tion with n > |, the amplitude of £ will be dominated 



by the C^2 mode 



/! ~ (ta -i>- 



3n 



(7) 



which is increasing. While for the expanding phase 
a ~ i n , Cfc is constant. This implies that for a cycle 
of cyclic universe during the contraction £fc is increased 
on large scale, up to the end of contracting phase in cor- 
responding cycle, while during the expansion it becomes 
constant. Thus the net result is that (k on large scale is 
amplified, which is inevitable. Here, n is constant only 
for the convenience of discussion. However, around the 
bounce, the case can not be described simply as Eq.(j6]), 
the qualitative result that the curvature perturbation on 
large scale is amplified is universal. 

The amplification for the perturbation modes, which 
are outside of Hubble scale in the j th and j th + 1 cycles 
all along, can be estimated as follows. In the j th cycle, 
after the bounce, (k is given by Eq.©, which will be un- 
changed up to the end of the j th cycle. Then the universe 
enters into the j th + 1 cycle, during the contraction of the 
j th + 1 cycle, (k w iU continue to increase. Thus at certain 
time during the contraction, in term of Eq.©, we have 
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since H ~ — t), which means that for these modes 

still staying outside of Hubble scale the amplitude of the 
spectrum will be amplified with same rate after each cy- 



cle. In the third line, 



tfj+l ~ (i 



+i 



ln( 



H. 



7+T 



which is 



J+i 



)ln 



H 3+i 
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FIG. 3: The evolutions of the comoving Hubble radius -ArT 
and the perturbation modes are plotted for the models with 
increasing cycles and equal cycles, respectively. 
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FIG. 4: The upper panel is the evolution of C,k in the model 
with increasing cycles, and the lower panel is that in the model 
with equal cycles. We can see that (k is increased cycle by 
cycle outside of Hubble scale, independent of the style of cyclic 
models. 
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is the efolding number for the primordial perturbation 
generated during the contraction of the j th + 1 cycle. 
Thus for the j th cycle, the amplification of the perturba- 
tion amplitude is ~ exp Af J . Thus with Eq. (jSJ) , we have 



t^ 1 - P+ 1 



1— 3nj+i 



n 



t j+i 

- 1 



Ce 
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where Q. is the perturbation generated during the first cy- 
cle. Eq. dTUl) indicates that the curvature perturbation is 
multipled through cycles, and thus the rate of the ampli- 
fication is quite rapid, which is general. We have assume 
that the linear perturbation approximation is satisfied 
all along. However, it might be broken after one or some 
cycles, which will bring a cutoff for j. 




FIG. 5: The evolution of $fc and S(j>k in the model with 
increasing cycles. We can see that 8<j>k is increased in con- 
tracting phase and constant in expanding phase, thus the net 
result is increased cycle by cycle. The corresponding mode is 

Kg 

In Fig. 4, the evolutions of the perturbation modes (k 
is plotted with the backgrounds in Fig.l by numerically 
solving Eqs.([T]),(I2]),(|3]). We can clearly see that is in- 
creasing during the contraction and is constant during 
the expansion, thus the net result is that the amplitude 
of perturbation is amplified cycle by cycle, which is con- 
sistent with Eq. (fT0|) . The metric perturbation $ and the 
perturbation 5(f) of background scalar field <j) are plotted 
in Figs. 5 and 6. They obviously also increase cycle by 
cycle. The results are consistent with above discussions, 
and also the numerical results in (2^1 , [H| > [HI for one cy- 
cle, and also 1411. 
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FIG. 6: The evolution of $fc and 8(j>k in the model with 
equal cycles. Their characters are similar with increasing 
cyclic models. The corresponding mode is 0.01fc s 



IV. THE EVOLUTION OF ISOCURVATURE 
PERTURBATIONS 

We will study the evolution of isocurvature perturba- 
tion. In principle, there may be a set of scalar fields 
together participating in the evolution of cyclic universe. 
The curvature perturbation is the perturbation along the 
field trajectory in fields space, while the isocurvature per- 
turbation is the perturbation orthogonal to the field tra- 
jectory [37[. The general case is slightly involved. How- 
ever, for simplicity, we assume that there is only a field 
X other than (f>, which has nothing to do with the back- 
ground evolution, and thus the background field <f> is re- 
sponsible for the curvature perturbation, while x omv 
contributes the isocurvature perturbation. 

The evolution of x field is given by 



X + 3H X + V x = 0. 



en; 



The simplest case is V(x) — 0, we have x ~ 1/a 3 - Thus 
if initially x 7^ 5 x w iU increase in the contracting phase, 
while x — > rapidly in the expanding phase. 
During the contraction phase ^ of each cycle, 



X 



dt 

^3 



(t B - t) 



l-3n 



(12) 



Thus x wm farther and farther deviate from its initial 
value in the model with equal cycles, while will eventually 
stop in certain value in the models with increasing cycle. 
Alternatively, the potential is 



V( X ) 



Mix 2 - 



(13) 
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During the contraction, the solution of Eq. (|TT|) is ap- 
proximately 



Mr, 



X 



3nM x (t B - 1) 



sm[M x (t B - t)}. (14) 



The field x will oscillate around its minimum x = and 
its amplitude will increase during the contraction of each 
cycle. 

The perturbation equation of x field is given by Eq. ([lj , 
6x k + 3H6 Xk + + V xx )5 X k = 4$ fe x - ^ k V x . (15) 



A. The evolution of Sx without the coupling to the 
metric perturbation 

We firstly neglect the coupling of x to the metric per- 
turbation <E>, which is right only when the effect of x on 
the background evolution is completely negligible. When 
the equation of Sx is 



k 2 

5xk + 3H5xk + ~^&Xk = 



V xx « H 2 



We define v = aSx, Eq. (IT6l) becomes 



[k 2 )v k = 0, 

a 



(16) 



(17) 



where the prime is the derivative for the conformal time 
■q. Initially, the perturbation is deeply inside its Hubble 



scale, fc 3> a!' I a, the solution is v k 



/2k 



e -ik V _ Dur _ 



ing the contraction, the perturbation is extended outside 
of the Hubble scale. When k 2 <C a" /a, the solution of 
Eq.tfTTf can be written as 



contraction in following cycle, its initial value is depen- 
dent on the ratio of the expansion of a to the contraction 
in previous cycle. In general, we have 

/ „• \ 3 



^Ee 
+ 1 



Ee 



n 



SXEe- (20) 



Thus if the cycle is equal, the initial value of Sxk is same 
for each cycle, however, if the cycle is increasing, it is 
decaying with the increasing of the number of cycles. 

The initial value of <5x^. +1 at the beginning time of the 
contraction determines the amplitude of the increasing 
of 5x^k +1 m the j th + 1 cycle. Thus after the bounce of 
the j th + 1 cycle, we generally have 



5 X k 



i+i 



!+ E Sx k +1 dt]Sxl (21) 



wheiedxi is the perturbation generated during the first 
cycle. Eq. flim) indicates that for the model with equal 
cycles, the perturbation modes will be amplified cycle 
by cycle and the amplification rate at each cycle is same, 
however, for the model with increasing cycles, since 8~x ~ 
1/a 3 is decaying with the increasing of the number of 
cycles, the integration J &x k +1 dt will be negligible after 
one or some cycles, which means that the perturbation 
Sxk will be amplified in the beginning and approach a 
constant with the increasing of the number of cycles. 

The evolutions of 5xk and in the models with equal 
cycles and increasing cycles are plotted in FigslT] and HI 
respectively. In FigJTl $Xk will increase cycle by cycle, 
and its increasing is slower than that of $k, which can 
be actually noticed by comparing Eqs. (fTU| and (|2Tj) . In 
FigE $Xk will only increase at the first cycle and then 
tends to constant in the following cycles. The numerical 
results are consistent with the analytic estimates. 



$Xk ^ C Sxl2 



dij 

la 



C s 



x,l 



(18) 



where the Cs x ,i mode is the constant mode and the Cs x ,2 
mode is the mode dependent on time. The spectrum of 
Sxk will be determined by which one of them is domi- 
nated. This solution is almostly same as ([5]) of curvature 
perturbation. The difference is that for £, the Cq^ mode 
is ~ J dr]/z 2 , which equals J dr)/a 2 only when lu is con- 
stant. However, around the bounce time, uj is rapidly 
changed. 

The increasing of the amplitude of Sxk during the con- 
traction is determined by the increasing mode 



dt 



1 



(tfl - t) 3n ~ 



(19) 
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FIG. 7: The solid line is the evolution of Sxk in the model 
with equal cycles, while the dashed line is the evolution of <!>fc. 



Thus Sxk ~ 1/a 3 . The matching condition around the 
bounce is that both Sxk and its time derivative are con- 
tinuous. Thus after the bounce 5x k will inherit the corre- 
sponding value before the bounce. However, since during 
the expansion 5~Xk ~ 1/a 3 is decreased. Thus before the 



B. The evolution of 8\ with the coupling to the 
metric perturbation 

During the contraction of each cycle, the energy of x 
field x 2 ~ 1/a 6 is rapidly increased. Its increasing is 
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FIG. 8: The solid line is the evolution of 8xt in the model 
with increasing cycles, while the dashed line is the evolution 
of 4> fe . 



generally faster than that of the background field 4>, thus 
it might be required to include the coupling of \ to the 
metric perturbation <£>. The metric perturbation is am- 
plified cycle by cycle, which is faster than that of free 6\, 
and thus will source the evolution of S\- We will check 
it. 

When V x = 0, that left in the right side of Eg. ([15]) is 
4$fcX, Eq. (fT5|) becomes 



3H6xk + -^6 X k = 4$fcX, 
Gr 



(22) 



In general, $ ~ H$> and \ — H\- What we care is the 
evolution of perturbation on large scale, i.e. k 2 <C a 2 H 2 . 
Thus Eq. (f22"jl becomes approximately 



and thus 



H 2 S X k ~ H^ kX , 



Sxk ~ 



(23) 



(24) 



In the model with equal cycles, a is periodically equal. 
In contracting phase of each cycle, we have 



1 



1 



(t B - ty- 



(25) 



where we have used a ~ (t B — t) n and H ~ 1/(*b — t). 
We can conclude that during the contraction with n > ^ , 
the increasing of the amplitude of Sxk is faster than that 
of the metric perturbation 

In the model with increasing cycles, x w i n become 
quite small after some cycles, since the maximal value 
of a in each cycle is increased with the increasing of the 
number of cycles. Thus after one or some cycles, will 
hardly have effect on Sxki therefore 6xu will approach a 
constant. 

Alternatively, for the potential (fT3)) . the equation of 8\ 

is 

S X k + 3H5 X k + (K + Ml)5 Xk = 4$ fc x - ^ k M 2 x (26) 



When M 2 <C H 2 , the result will be same with that of 
V = 0. While for M 2 » ff 2 , we have 

S X k - X$fc, (27) 

where x is given by Eq. (fT4"]) . Thus the solution (|27|) is 
approximately 

5 * k ~ m !f P A MM x (t B - t)]$ k . (28) 
M x {t B - t) 

Thus the perturbation will oscillate and its amplitude 
will be amplified faster than However, the increasing 
rate of the amplitude dose not decay with the increasing 
of the number of cycles. 
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FIG. 9: The background is that of the model with equal 
cycles. The solid line is the evolution of 4> fc and the dashed 
line and the dotted line are the evolution of Sxk with V(x) ~ 
and V(x) = iM^x 2 , respectively. 
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FIG. 10: The background is that of the model with increasing 
cycles. The solid line is the evolution of and the dashed 
line and the dotted line are the evolution of Sxk with V(x) = 
and V(x) = |M^x 2 , respectively. 



The numerical results for the models with equal cycles 
and increasing cycle, respectively, are plotted in Figsj9] 
andHUl We can see that the perturbation $Xk is amplified 
faster than cycle by cycle. The oscillating of x leads 
the oscillating of Sxk- These results are consistent with 
analytical discussions. 



V. DISCUSSION 

The cosmological cyclic scenario offered an alternative 
insight into the origin of observable universe. In general, 
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it is thought that the background evolution of cyclic uni- 
verse is globally homogeneous cycle by cycle. However, 
the amplitude of curvature perturbation on large scale 
will be amplified cycle by cycle. 

In certain sense, the amplification of curvature per- 
turbation might eventually destroy the homogeneity of 
background, which might lead to the ultimate end of cy- 
cles of global universe. However, it can be argued (l5lllq] 
that for the model with increasing cycles, the global uni- 
verse will possibly evolve into a fissiparous multiverses 
after one or some cycles, in which the cycles will con- 
tinue only at corresponding local regions, inside of which 
the background is homogeneous. 

The amplitude of the isocurvature perturbation on 
large scale may be also amplified cycle by cycle. How- 
ever, the rate of its amplification is slower than that of 



curvature perturbation, and further for the model with 
increasing cycles, with the increasing of the number of cy- 
cles the amplification of its amplitude will be suppressed, 
and after one or some cycles the resulting amplitude will 
approach a constant. However, if the coupling to the met- 
ric perturbation is included, the result will be altered. 
This might have interesting implication to the applica- 
tion of the curvaton mechanism [i^] , [H[ , [3 , [H| , called 
the bounce curvaton scenario |46| . 
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